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ABSTRACT 

We present an analysis of XMM-Newton X-ray data of WR 30a (WO+O), a close 
massive binary that harbours an oxygen-rich Wolf-Rayet star. Its spectrum is charac¬ 
terized by the presence of two well-separated broad peaks, or ‘bumps’, one peaking at 
energies between 1 and 2 keV and the other between 5 and 7 keV. A two-component 
model is required to match the observed spectrum. The higher energy spectral peak 
is considerably more absorbed and dominates the X-ray luminosity. For the currently 
accepted distance of 7.77 kpc, the X-ray luminosity of WR 30a is Lx > 10 34 erg s _1 , 
making it one of the most X-ray luminous WR+O binary amongst those in the Galaxy 
with orbital periods less than ~ 20 d. The X-ray spectrum can be acceptably fitted 
using either thermal or nonthermal models, so the X-ray production mechanism is yet 
unclear. 
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1 INTRODUCTION 

Wolf-Rayet (WR) stars are massive stars that have pow¬ 
erful winds and are losing mass at high rates (M ~ 10~ 5 
Mq yr -1 ; V w i n d = 1000 — 5000 kms -1 ). Based on their op¬ 
tical spectra, they are divided into three subtypes: nitrogen- 
rich (WN), carbon-rich (WC), and oxygen-rich (WO). Their 
progenitors are massive stars with initial masses >25 Mq 
and most WR stars are thou ght to end their lives as su¬ 
pernovae (see ICrowtheil 120071 for a review on the physical 
properties of WRs). The observed binary fraction in Galac¬ 
tic Wolf-Rayet stars is relativ ely high with about 4 0% being 
members of WR+O systems (Ivan der Huchtll200ll 'l. 

WR stars were di scovered to be X-ra y sources by the 
Einstein Observatory (1 Seward et all Il979k The first sys¬ 
tematic survey of WRs showed that WR+O b inaries ar e 
the brightest X-ray sources amongst them dPollockl 1 198715 . 
Their enhanced emission likely originates from the interac- 
tion region of the winds of the massive binary components 
d Prilutskii fe Usovlll976l : ICherepashchuk|[l976l l . 

After the launch of the modern X-ray observatories 
Chandra and XMM-Newton , the number of WR stars with 
good quality X-ray spectra has increased considerably. Thus, 
some similarities and differences in the characteristics of 
their X-ray emission are now established. 

Most WR stars are classified as either WC or WN sub- 
types and analysis of undispersed spectra of presumably sin¬ 
gle objects have revealed clear differences between WC and 
WN stars. All observations of putati vely single WC stars 
so far have resulted in non-detections (lOskinova et al.ll2003l : 
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ISkinner et alJlioObl l. In contrast, WN stars without known 
companions are detected in X-rays and their spectra typi¬ 
cally reveal an ad mixture of cool (kT < 1 keV) and hot (kT 
> 2 keV) plasma (ISkinner et al.ll2010l . 2012). WO stars are 
much rarer. The only pointed X-ray observations of a pre¬ 
sumably single Galactic WO star to date, WR 142, res ulted 
in a detection dOskinova et al.l 120091 : ISokal et al.1 l20ld l . In¬ 
terestingly, the Chandra detection of WR 142 showed an 
extremely hard and heavily-absorbed X-ray spectrum which 
coul d be fitt e d us ing either thermal or nonthermal models 
( Sokal et al.l l201Clh . 

On the other hand, the X-ray properties of massive 
WR+O binaries do not seem to have any strong dependence 
on the subtype of the WR component. Analyses of the X- 
ray emission from both close and wide WR+O binaries give 
support to the idea that a substantial fraction of their X- 
ray emission arises in colliding stellar wind (CSW) shocks, 
resulting from the inter action of the massive winds of the 


binary components (e.g. 

Skinner et al. 2001 

: Raassen et al. 

2003); 

Scliild et al. 2004]; 

Pollock et al.l 2005|; 

Zhekov & Park 

2010a 

. b: Zhekov et al. 

2Oil], 2014; Zhekov 

l2012h. But no 


WR+O binary with an oxygen-rich WR component has 
heretofore been detected in X-rays. 


In this paper, we report the X-ray detection of the close 
WO+O binary system WR 30a. The paper is organized as 
follows. We summarize information on WR 30a in Section [2] 
In Section [3] we review the X-ray observations and data. 
In Section [4] we present results from analysis of the X-ray 
properties of WR 30a. In Section [5] we discuss our results, 
and we present our conclusions in Section [6] 
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Figure 1. Top row. The total-field high-energy (>10 keV) EPIC light curves of WR 30a. Middle row. The total field high-energy EPIC 
light curves of WR 30a after being filtered to remove high X-ray background. Bottom row. The filtered EPIC images of WR 30a in the 
0.2-10 keV energy band. The source spectrum was extracted from within the central circle (radius of 30"). The background spectrum 
was extracted from the surrounding annulus (outer radius of 60 ,/ ). The circled plus sign gives the optical position of WR 30a (SIMBAD). 


2 THE WOLF-RAYET STAR WR 30a 

WR 30a (V574 Car; [MS70]) is o ne o f four known oxygen- 
rich WR stars in the Galaxy (Ivan der Huchtl l200ll . 2006; 
also Galactic Wolf-Rayet Catalogu^H). It is a W04+05- 
5.5 b inary system at a distance of 7.77 kpc (Ivan der Huchtl 
l200ll) . Detailed photometric and spectroscopic studies pro¬ 
posed an orbital period of 4.619 d and classified the 05 
component as either _a main-sequence or, more likely, a 
giant star dGosset et al.l l200lh . The optical extinction to- 
ward WR 30a is A„ = 4.26 mag llvan der Huchtl l200ll : 
A„ = 1.11 Ay) implying a foreground column density 
of Nh = (6.33 — 8.52) x 10 21 cm -2 . The range corre¬ 
sponds to the conversion that is used: N h = (1.6 — 1 .7) x 
10 21 Ay cm -2 (lYuong et al.llio lid. ICetrna n et al.l l2005l ): and 
N h = 2.22 x 10 21 Ay cm -2 ( Gorensteinl 19750 . Similarly 
to other WO stars, WR 30a has a very fast stellar wind: 
V w i„d = 4500 kms~ 1 dKingsburgh et al.ll l 9951 ). 

No radio data o n WR 30a were found in the literature. 
IPollock et al „ d 199.il l reported an X-ray non-detection in the 


1 http://pacrowther.staff.shef.ac.uk/WRcat/index.php 


ROSAT PSPC Survey on Wolf-Rayet stars, likely due to 
the very short exposure time (260 s). 


3 OBSERVATIONS AND DATA REDUCTION 

We searched the Chandra and XMM-Newton archives for 
data on WR 30a and found one X-ray data set. WR 30a was 
observed with XMM-Newton on 2010 January 25 (Observa¬ 
tion ID 0606330101) with a nominal exposure of ~ 60 ks. 

Our study is based on the data from the European Photon 
Imaging Camera (EPIC), which has one pn and two nearly- 
identical MOS detector^. We used the XMM-Newton SA10 
14.0.0 data analysis software for data reduction, following 
the procedures summarized below. 

First, the SAS pipeline processing scripts emproc and 
epproc were executed to incorporate the most recent calibra¬ 
tion files (as of 2015 January 27). The data were checked for 
high X-ray background following the instructions in the SAS 

2 See § 3.3 in the XMM-Newton Users Hand¬ 
book, Issue 2.12, 2014 (ESA: XMM-Newton SOC) 

http://xmm.esac.esa.int/external/xmm_user_support/documentation/uhb 
a Science Analysis Software, http://xmni.esac.esa.int/sas 
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documentatior@, namely, by constructing total-field light 
curves at high energies: E > 10 keV for M0S1,2 and E in the 
10 - 12 keV range for pn. As seen from Fig. |T] a high (and 
flaring) X-ray background was registered in the EPIC ex¬ 
posure, with the pn data being most affected. Nevertheless, 
we applied a background filtering procedure that allowed us 
to make use of some of the pn exposure and a substantial 
fraction of the MOS exposures. 

We made use of the XMM-Newton Extended Source 
Analysis Software (xmm-esas) packag43, which is now in¬ 
corporated in SAS. We created filtered event files using the 
ESAS commands mos-filter and pn-filter which minimize the 
contamination of the soft proton flaring in the EPIC data in 
a robust manner (for details see the xmm-esas CookboolQ). 
As a check, we constructed the total-field filtered light curves 
at high energies (Fig. [lj which provided confidence that the 
background had been reduced to an accepable level. 

We then proceeded with the source and background 
spectral extraction from the pn and MOSl,2 filtered event 
files using the extraction regions shown in Fig. [l] The SAS 
procedures rmfgen and arfgen were used to generate the cor¬ 
responding response matrix files and ancillary response files 
for each spectrum. The MOS spectrum used in our anal¬ 
ysis is the sum of the spectra from the two MOS detec¬ 
tors. The extracted spectra (0.2 - 10 keV) of WR 30a had 
~ 403 source counts (~ 710 source+background counts) in 
the 24.1-ks pn effective exposure and ~ 350 source counts 
(~ 584 source+background counts) in the 37.2-ks MOS ef¬ 
fective exposure (36.0 ks in MOS1 and 38.4 ks in MOS2). 

Since the pn exposure was more heavily affected by high 
X-ray background, our analysis emphasized the MOS spec¬ 
trum while the pn spectrum was used to check the con¬ 
sistency of the results. For the spect ral analysis, we used 
version 12.8.2 of xspec l|Arnau dl ll996h . 

Also, we constructed the pn and MOSl,2 background- 
subtracted light curves of WR 30a. Since the deleted high- 
background time intervals were unevenly-spaced, the good 
time intervals (GTIs) in the filtered light curves were as well. 
There were 27 GTI segments as short as 120 s and as long as 
6,540 s in MOS1; 38 GTI segments as short as 120 s and as 
long as 4,315 s in MOS; and 26 GTI segments as short as 15 
s and as long as 7,080 s in pn. On a timescale less than 24-37 
ks, the X-ray light curves were statistically consistent with 
a constant count rate. A constant count rate model gives a 
goodness-of-fit ^ 0.98 using \ 2 fitting. 


4 RESULTS 

As apparent from Fig. [2] the X-ray spectra are nearly 
featureless and are characterized by two broad peaks or 
‘bumps’, one at lower energies (peaking between 1 and 2 
keV) and another at higher energies (maximum emission 
between 5 and 7 keV). Other than a possible contribution 
from the Fe K complex near 6.67 keV, there are no bright 
spectral lines that would clearly indicate thermal emission. 


We thus explored both thermal (optically-thin plasma) and 
non-thermal (power-law) models in XSPEC. 

For the ther mal models, we adopted ty pical WO star 
abundances from Ivan der Hucht et al.1 (j 19866 . The emission 
plasma components were subject to common X-ray absorp¬ 
tion (model wabs in xspec) representative of the interstel¬ 
lar medium (ISM). Some additional absorption, as expected 
from the winds, was also required to acceptably fit the total 
spectrum. Typical WO star abundances were adopted for 
the wind absorption which was modeled using phabs model 
in xspec. Due to the rather low number of X-ray counts and 
absence of strong lines, all the abundances were kept fixed 
at their adopted values during fitting. 

The one-component model subject only to ISM X-ray 
absorption did not provide a good fit to the observed MOS 
spectrum of WR 30a. The same was true for the two- 
component model when using the same ISM absorption for 
each component. 

The two-component model with individual X-ray ab¬ 
sorptions for each emission component gave a very good 
fit to the MOS spectrum. However, the best-fit plasma- 
temperature (or power-law index) values of the different 
emission components overlapped at the lcr level. We thus 
adopted a common value for these parameters in the final 
model fits, which means that the spectral shape of the emis¬ 
sion components was the same but the two components dif¬ 
fered in their respective amounts of thermal plasma (or non- 
thermal emission) and in their X-ray absorptions. When ex¬ 
ploring the cases with wind absorption, we kept the common 
ISM X-ray absorption fixed to the values that correspond to 
the optical extinction of WR 30a (Section 0. As an addi¬ 
tional check of our global modelling, we fitted the MOS and 
pn spectra of WR 30a simultaneously by making use of the 
same models described above. 

Table [T] and Fig. [5] present the corresponding results 
from the fits to the XMM-Newton spectra of WR 30a. The 
fit results can be briefly summarized as follows. Both ther¬ 
mal and non-thermal models are successful in reproducing 
the X-ray spectrum. In the case of thermal models, a rela¬ 
tively hot plasma (kT > 2 keV) is required to be present in 
the X-ray emitting region. The first emission component in 
the models is subject to X-ray absorption that is in general 
consistent with the optical extinction toward WR 30a, while 
the second emission component is more heavily absorbed. If 
we adopt a common ISM absorption for both emission com¬ 
ponents, then only the second component requires some ad¬ 
ditional wind absorption (see Models B and D in Tabic |T]). 
For illustration of the results from the model fits with a com¬ 
mon ISM absorption, we show only the case with the upper 
limit for the range of values that correspond to the optical 
extinction (Section^. We note that if the lower limit for the 
ISM absorption was used, then the derived model parame¬ 
ters were within their respective lcr confidence intervals. 


4 see § 4.4.4 in Users Guide to the XMM-Newton Science The X-ray luminosity of WR 30a , Lx > 10 34 erg s Sis 

Analysis System, Issue 11.0, 2014 (ESA: XMM-Newton SOC); ^ite hi S h ™ hen compared to other WR stars in the Galaxy. 

http://xmm.esac.esa.int/external/xmiiLuser.support/documentationyilcisiHiy^/lf&iJ/'ily absorbed component responsible for the 

5 http://xmm.esac.esa.int/sas/current/doc/esas observed emission in the 5-7 keV range dominates the total 

6 ftp://xmm.esac.esa.int/pub/xmm-esas/xmm-esas.pdf luminosity in the 0.5-10 keV energy band. 
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Table 1. Global Spectral Model Results 
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Note. Results from fits to the EPIC spectra of WR 30a. The labels MOS and MOS+pn denote which spectra were fitted. The 
XSPEC models are: wabs * apec + wabs * apec (Model A); wabs * ( phabs * apec + phabs * apec) (Model B); wabs * powerlaw-\- 
wabs * powerlaw (Model C); wabs * (phabs * powerlaw + phabs * powerlaw) (Model D). Note that both model components have the 
same spectral shape (see text for details). Tabulated quantities are the neutral hydrogen absorption column density (N#), helium 
absorption column density (Njj e ) for the ‘wind’ components, plasma temperature (kT), emission measure (EM = f n e n# e dV), photon 
power-law index (Tpoty), power-law model normalization (norm pott; ), the observed X-ray flux (Fx) in the 0.5 - 10 keV range followed in 
parentheses by the unabsorbed value and the X-ray luminosity (Lx)- Fx l and Lx.i denote the values for the first model component, 
peaking at lower energies. The abudnaces for the thermal models are those typical for the WO stars (by number: H = 0.00, He = 0.266, 
C = 0.213, N = 0.00, O = 0.497, Ne = 9.69 x 10“ 3 , Mg = 1.26 x 10 -2 , Si = 6.73 x 10 -4 , S = 1.47 x 10 -4 , Fe = 3.69 x 10 -4 ; 

Ivan der Hucht et al.|[l986r) . The values for the emission measure and the X-ray luminosity are for a reference distance of d = 7.77 kpc. 

Errors are the lcr values from the fits. 
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Figure 2. The background-subtracted EPIC spectra of WR 30a overlaid with some best-fitting models (Table [T]). In the case of 
simultaneous fit to the MOS and pn spectra, the latter is plotted in red. For presentation purposes, the spectra were slightly re-binned 
with respect to the original binning of 20 counts per bin used in the fits. The model components are shown in blue and green colour. 
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5 DISCUSSION 


Some of the results from the X-ray spectral fits may be af¬ 
fected by the fact that there are relatively few counts in 
the XMM-Newton spectra (Section^. Nevertheless, the fol¬ 
lowing physical picture emerges. There is likely some distri¬ 
bution of thermal X-ray plasma (or non-thermal emission) 
with similar plasma characteristics and this X-ray emission 
region is only partially absorbed. Alternatively, there may 
be two different X-ray sources with about the same plasma 
characteristics but one of them is more heavily absorbed. If 
the absorption is due to the stellar wind, that source must be 
located closer to the WO star while the other less-absorbed 
source is further from the star and is not subject to signifi¬ 
cant wind absorption. 

If the X-rays from WR 30a are from CSW shocks in the 
WR+O binary, we expect a temperature-stratified interac¬ 
tion region. Thus, the derived values of the plasma tempera¬ 
ture in our fits represent a mean temperature for the shocked 
plasma. For the case of typical WO abundances, the post¬ 
shock plasma temperature is kT = 3.36 Uiooo keV, where 
Viooo is the shock velocity in units of 1000 kms -1 . Given the 
stellar wind velocity of 4500 km s -1 (Section[5J , high temper¬ 
atures at kT > 2 keV are realistic for shocks in WO plasmas. 
Very high shock temperatures could be achieved even if the 
winds have not reached their terminal velocities before col¬ 
liding, which could be the case in this close WO+O binary 
(orbital period of 4.619 d; Section [2]). 

Even though the high plasma temperatures (kT > 2 
keV) determined from spectral fits could be taken as evi¬ 
dence for colliding stellar winds, it is worth keeping in mind 
that such high temperatures have also been dete cted in pre¬ 
sumably single WN stars (|Skinner et ah 2010|, 2012 ) and 
in very close binaries like WR 46 (see Gosset - j^d. 201l|; 
Zhekov] 120121 1 and WR 155 (CQ Cep; see ISkinner et all 


2015l l. So. CSWs may not be the only mechanism producing 


very hot plasma in WR stars. 

However, as already noted, the spectral fits of 
WR 30a cannot clearly distinguish between thermal and 
nonthermal models. This ambiguity is due in part to the low 
number of spectral counts and higher signal-to-noise spec¬ 
tra capable of determining whether weak emission lines are 
present will be needed to discriminate between thermal and 
nonthermal emission. A similar ambiguity was found for the 
presumably single WO star WR 142 wh ose X-ray emissio n 
is even fainter and shows no strong lines JSokal et al.ll2010f) . 
WR 30a and WR 142 also reveal other similarities, e.g. both 
have a hard heavily-absorbed spectral component. A notable 
difference is that the X-ray luminosity of WR 142 is less than 
1% of that of WR 30a. The binary nature of WR 30a may be 
largely responsible for this, but distance uncertainties could 
also contribute. 


The X-ray luminosity of WR 30a (Lx > 10 34 erg s _1 , 
Table[lJ may be the highest amongst the close WR+O bina¬ 
ries (orbital pe riods less than ~ 20 d) in the Galaxy studied 
so far (e.g., see lZhekovll2012l : ISkinner et al.ll2015l : lNaze et al.l 
2008), provided the adopted distance of 7.77 kpc is not sig¬ 
nificantly overestimated. The high X-ray luminosity must be 
accommodated by any plausible theoretical emission mod¬ 
els. 


We note that other close WR binari es with very high X - 
ray luminosities may exist in the Galaxy. I Gagne et al.l (120121 1 


have compiled a list of possi- ble detections of colliding wind 
binaries with Chandra , XMM-Newton and ROSAT. High 
X-ray luminosity candidates are WR28, WR29, WR43a, 
WR43c and WRlOlk. However, caution is needed when con¬ 
sidering the available data on these objects. For example, 
WR43a and WR43c are located in the core of the Galac¬ 
tic starburst region NGC 3603 and thus cannot be spatially 
resolved in X-rays (that is their X-ray spectra are contam¬ 
inated by the near-by sources). Similarly, WRlOlk cannot 
be resolved since it is in the centre of the Galaxy only 1.8" 
from Sgr A (SIMBAD). The data on WR28 and WR29 are 
from ROSAT observations and only upper limits can be de¬ 
rived on their observed fluxes. Nevertheless, these close WR 
binaries which may have extremely high X-ray luminosities 
deserve more attention and in-depth X-ray studies. 

Finally, a distinguishing feature of WR 30a is its un¬ 
usual spectral shape with two broad separated peaks. After 
examining the X-ray spectra of other WR stars we were able 
to identify another object with a similarly shaped spectrum, 
the close WR+O binary WR 79. 


5.1 Comparison with WR 79 

WR 79 is a WC+O binar y with an o rb ital pe riod of 8.28 d at 
a distance of 1.99 kpc (Ivan der Hu cht 2001). Figure [ 3 ] com¬ 
pares its Chandra X-ray spectrum with that of WR 30a and 
WR 142. We note that WR 79 was also observed by XMM- 
Newton but its spectrum is contaminated by two nearby 
sources within 7” - 8" that are seen in the higher spatial 
resolution C handra image (for details, see Appendix A in 
IZhekovll20lll . 

Two broad peaks are present in both WR 79 and 
WR 30a, but only the higher energy peak is present in the 
apparently single WO star WR 142. No prominent spectral 
lines are seen with the possible exception of Fe emission 
in the 6-7 keV range in WR 79 and a hint of the same in 
WR 30a. But faint emission lines could be missed since the 
number of counts in all three spectr a is quite lim ited, with 
only ~ 380 source counts in WR 79 d Zlieko vl 1201 2l l and ~46 
counts in WR 142 (ISokal et al.ll2Q10l l. 

The X-ray spectrum of WR 79 was well matched by a 
two-component optically-thin plasma model with individua l 
wind absorptions (see Table 3 and Fig. 1 in IZhekovll2012l 'l. 
As the global fits showed, the second model component due 
to hot plasma with kT > 2 keV is subject to much higher 
wind absorption, which is also the case with WR 30a. And 
the more absorbed component, whose emission produces the 
broad peak between 4 and 7 keV in the spectrum, dominates 
the total luminosity in the 0.5-10 keV energy range. We did 
this estimate by making use of exactly th e same s pectra l 
data (no re-processing) and models as in IZhekovI (1201214 . 
Despite these similarities, there is an interesting difference: 
the X-ray luminosity of WR 79 is less than WR 30a by 
approximately one order of magnitude. 


6 CONCLUSIONS 

In this work, we presented and analyzed archival XMM- 
Newton data of WR 30a which provide the first X-ray spec¬ 
tra of this WR+O binary that harbours a rare oxygen-rich 
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Figure 3. The background-subtracted spectra of 

WR 30a (XMM- Newton EPIC -MOS), WR 79 (Chandra ACIS-S, 
Obsld 5372, see IZhekovI 1201 and WR 142 (Chandra ACIS-I, 
Obsld 9914; based on our own spectral extraction using version 
4.7 of Chandra Interactive Analysis of Observations, CIAO; 
http://cxc.harvard.edu/ciao/). 

Wolf-Rayet star. The main results and conclusions are as 
follows. 

(i) The X-ray spectrum of WR 30a shows no strong 
spectral lines, although some Fe emission between 6-7 keV 
may be present and other weak lines may have been missed 
due to the low number of source counts in the spectrum. 
The spectrum can be acceptably fitted by either thermal or 
nonthermal emission models. 

(ii) The most distinguishing feature in the WR 30a spec¬ 
trum is the presence of two broad but well-separated peaks, 
one peaking at energies between 1 and 2 keV and the other 
between 5 and 7 keV. A two-component model is required to 
match the observed spectrum. The spectral component seen 
in the higher energy peak is considerably more absorbed and 
dominates the intrinsic X-ray luminosity in the 0.5-10 keV 
range. 

(iii) For the currently accepted distance of 7.77 kpc, the 
X-ray luminosity of WR 30a is Lx > 10 34 erg s _1 , making it 
the most X-ray luminous object amongst the close WR+O 
binaries (orbital periods less than ~ 20 d) in the Galaxy 
studied so far. 

(iv) On the basis of the double-peaked spectral shape 
and other X-ray characteristics, we have identified two 
WR binaries that are quite similar: WR 30a and WR 79 
(WC+O). 

(v) Future observations of both single and binary WO 
stars are needed to obtain higher signal-to-noise X-ray spec¬ 
tra capable of determining if faint emission lines are present 
and thereby distinguishing between thermal and nonthermal 
emission processes. 
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